Approximate bound state solutions of the deformed 
Woods-Saxon potential using asymptotic iteration 

method 

Babatunde J. Falaye(^ 

Theoretical Physics Section, Department of Physics 
University of Ilorin, P. M. B. 1515, Ilorin, Nigeria. 



Majid Hamzavi[^ 

Department of Basic Sciences, Shahrood Branch, 
Islamic Azad University, Shahrood, Iran 



Sameer M. Ikhdairlf] 

Physics Department, Near East University, 
922022 Nicosia, North Cyprus, Mersin 10, Turkey 



Abstract 

By using the Pekeris approximation, the Schrodinger equation is approximately solved for the 
nuclear deformed Woods-Saxon potential within the framework of the asymptotic iteration 
method. The energy levels are worked out and the corresponding normalized eigenfunctions 
are obtained in terms of hypergeometric function. 
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1 Introduction 

The deformed Woods-Saxon (dWS) potential is a short range potential and widely used in nu- 
clear, particle, atomic, condensed matter and chemical physics [1-7]. This potential is reasonable 
for nuclear shell models and used to represent the distribution of nuclear densities. The dWS 
and spin-orbit interaction are important and applicable to deformed nuclei [8j and to strongly 
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deformed nuclides [9j. The dWS potential parameterization at large deformations for plutonium 
237,239,241 odd isotopes was analyzed [10]. The structure of single-particle states in the second 
minima of 237,239,241 ^as been calculated with an exactly dWS potential. The Nuclear shape 
was parameterized. 

The parameterization of the spin-orbit part of the potential was obtained in the region corre- 
sponding to large deformations (second minima) depending only on the nuclear surface area. The 
spin-orbit interaction of a particle in a non-central self consistent field of the WS type potential 
was investigated for light nuclei and the scheme of single-particle states has been found for mass 
number Aq = 10 and 25 [8j . Two parameters of the spin-orbit part of the dWS potential, namely 
the strength parameter and radius parameter were adjusted to reproduce the spins for the values 
of the nuclear deformation parameters |llj . 

Badalov et al. investigated Woods-Saxon potential in the framework of Schrodinger and 
Klein-Gordon equations by means of Nikiforov-Uvarov method \12\ I13j . In our recent works 
\14:\ [T5] , we have studied the relativistic Duffin-Kemmer-Petiau and the Dirac equation for dWS 
potential. In addition, we have also obtained the bound state solutions of the PT — /non — PT- 
symmetric and non-Hermitian modified Woods-Saxon potential with the real and complex- valued 
energy levels |16j . 

In this paper we will study Schrodinger equation with dWS potential for any arbitrary orbital 
quantum number We obtain the analytical expressions for the energy levels and wave functions 
in closed form. Therefore, this work is arranged as follows: In section 2, a brief introduction to 
asymptotic iteration method (AIM) is given. In section 3, the Schrodinger equation with dWS 
potential is solved. Finally, results and conclusions are presented in section 4. 

2 The Asymptotic Iteration Method 

We briefly outline the AIM here; the details can be found in references |17l I18j . 
2.1 Energy eigenvalues 

AIM is proposed to solve the homogenous linear second-order differential equation of the form 

y'nix) = Xo{x)y'n{x) + So{x)yn{x), (1) 

where Xoix) ^ and the prime denotes the derivative with respect to x, the extral parameter n is 
thought as a radial quantum number. The variables, So{x) and Xo{x) are sufficiently differentiable. 
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To find a general solution to this equation, we differentiate equation ([T]) with respect to x as 

y'nix) = Ai(3;)y^(x) + si{x)yn{x), (2) 

where 

Ai(x) = \'^{x) + So{x) + Xl{x), 

si{x) = s'^{x) + So{x)\o{x), (3) 
and the second derivative of equation ([T]) is obtained as 

y'nix) = X2{x)y'n{x) + S2(a;)y„(x), (4) 

where 



A2(x) = X'i{x) + si{x) + Xo{x)Xi{x), 
S2ix) = s'i{x) + So{x)Xi{x). 



(5) 



Equation Q can be iterated up to {k + l)th and {k + 2)th derivatives, k = 1,2, 3... Therefore we 
have 



where 



Vn^^Hx) = Xk-i{x)y'^{x) + Sk-i{x)yn{x) 
y'^~^^\x) = Xk{x)y'n{x) + Sk{x)yn{x), 



Xk{x) = X'k_i{x) + Sk-i{x) + Xo{x)Xk-i{x), 

Sk{x) = Sk^i{x) + So{x)Xk-l{x). 



(6) 



(7) 



From the ratio of the (k+2)th and (k+l)th derivatives, we obtain 



dx 



-In 



(fc+2). 



Xk{x) 


y'n{x) + y^^\]-^yn{x)] 




Xk-i{x) 


[y'ni^) + x\_\[x)yn{x)\ 



(8) 



if A; > 0, for sufficiently large k, we obtain a values from 

Sk{x) Sk-l{x) 



with quantization condition 

5k{x) 



Afc(x) Xk-i{x) 



Afc(x) Sk{x) 
Afc_i(x) Sk-i{x) 



a{x). 



, k = 1,2,3.. 



(9) 



(10) 
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Then equation (p| reduces to 



x] 



Afc(x) 



which yields the general solution of Eq. ([T]) 

Unix) = exp ^— J a{x')dx'^ C2 + Ci J exp ( J [Xo{x") + 2a{x")] dx" j dx' 



(11) 



(12) 



For a given potential, the idea is to convert the radial Schrodinger equation to the form of 
equation ([T]). Then Xo{x) and So{x) are determine and Sk{x) and Xk{x) parameters are calculated 
by the recurrence relations given by equation ([T]). The energy eigenvalues are then obtained by 



the condition given by equation (10) if the problem is exactly solvable. 



2.2 Energy eigenfunction 

Suppose we wish to solve the radial Schrodinger equation for which the homogenous linear second- 
order differential equation takes the following general form 



y"{x) 



tx 



N+l 



m + 1 



X 



y'{x) 



Wx 



N 



1 - 6x^+2 ■ 



\l-bx^+^ 

The exact solution y„,(x) can be expressed as [18] 

Vnix) = i-irC2{N + 2)"(a)„2i^i(-n, p + n;a; 



where the following notations has been used 

2m + TV + 3 



r((j + n) 



3 Any ^-state solutions 



A^ + 2 



and p 



{2m + l)b + 2t 
{N + 2)b 



(13) 



(14) 



(15) 



The deformed Woods-Saxon potential we investigate in this study is defined as [HI [151 [E] 

V{r) = -- ^^TTTPV, i? = ro^y^ Vo = {40. 5 + 0.13 Ao)MeV, R»a, g > 0, (16) 

1 -hgexpC^) 

where Vq is the depth of potential, g is a real parameter which determines the shape (deformation) 
of the potential, a is the diffuseness of the nuclear surface, R is the width of the potential, j\.q is 
the atomic mass number of target nucleus and Tq is radius parameter. By inserting this potential 
into the Schrodinger equation [T71 [18] as 



'2/x 



I d . d 



1 



1 



d / ^ d 

^2 Qj. J.2 g^j^ do \ 89 J sin^ 9 dcj)^ 



(17) 



and setting the wave functions ^'nfm(^) = -Rn£(^)^m(^) </')^ ^) we obtain the radial part of the 
equation as 

r ^2 9/, / \ i(£ + l)' 

'^'^ > Rni{r) = 0. (18) 



^2 2n f 



l + gexp(^)^ 



Because of the total angular momentum centrifugal term, equation (18) cannot be solved ana- 



lytically for i ^ 0. Therefore, we shall use the Pekeris approximation in order to deal with this 
centrifugal term and we may express it as follows [201 EI] 
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Ucent. (r) 



1 



+ 



with X = r — R. In addition, we may also approximately express it in the following way 



(19) 



U 



1 



1 



+ 



Do 



(20) 



l + gexp(i/x) ' (1 + gexp (i/x))^ 
where v = 1/a. After expanding (20) in terms of x, x^, x^,- • • and next, comparing with equation 



(19), we obtain expansion coefficients Do, Di and D2 as follows: 



Do 
Di 
D2 



i-il±il!(i 



3 



i^Rq^ 
uRq^ 



-1 + 3g 



6(1 + g) 
vR 



(21) 



l-q , 3(l+(?) 



+ 



uR 



Now, inserting the approximation expression (20) into equation (18) and changing the variables 



z through the mapping function z = (1 + gexp {vx)) ^, equation (18) turns to 



d'Rneiz) ^ l-2z dRnejz) ^ 



z{l — z) dz 
Do£(i + 1 



[uz{i - z)r 



2fiVo ^ Die{£ + 1) 



D2zHii + 1) 



Rne{z) 



0.(22) 



R^ J \ i?2 J J12 

Before applying the AIM to this problem, we have to obtain the asymptotic wave functions 



and then transform equation ( 22 ) into a suitable form of the AIM. This can be achieved by the 



analysis of the asymptotic behaviours at the origin and at infinity. As a result the boundary 
conditions of the wave functions Rni{z) are taken as follows: 



Rne{z) —7- when z — )• 1, 
Rne{z) — ?• when z — )• 0, 



(23) 



thus, one can write the wave functions for this problem as 



Rni{z) = z^{l-zr'F^,{z), 



(24) 



where we have introduced parameters a and 7 defined by 



a 



7 



i?2 



n 



+ 1) 



{Do + D1+D2 



(25) 



for simphcity. By substituting equation ( 24 ) into equation ( 22 ) , we have the second-order homo- 



geneous differential equation of the form: 
■(2a + l)-2z(a + 7 + l 



Keiz) + 



zil-z) 



Keiz) 



a + 7 + a + 7 + 2a7 



z{l-z) 



Fnl{z) = 0, 

(26) 

which is now suitable to an AIM solutions. By comparing this equation with equation ([T]), we 
can write the \o{z) and So{z) values and consequently; by means of equation ([T]), we may derive 
the Afc(z) and Sk{z) as follows: 



\o{z) 
So{z) 

Xliz) 



2z(a + 7 + l) - (2a + 1) 



z{l-z) 



+ 7^ + a + 7 + 2q7 



z(l-z) 



2 + z (2 - 37 + 7 - ^TiTF 



2 D2e{e+ 



11) _ ^2 _ 3^2 



D2£(m) 



z2(l -z)2 

[a{z - l)(z(67 - 7) - 6)] + 0^(4 -7z + 3z^) 



(27) 



si{z) 
. . etc. 



(a + 7)2 + a + 7 ■ 



bz-il + ail-z))] 



z2(l -z)5 



The substitution of the above equations into equation (|10|), we obtain the first 6 values as 
(a + 7)2 + (a + 7) 



6o{z) 



1/2 i?2 



a(3 + 27) + 7(3 + 7) -2 + a 



2 _ D2£{m) 



z2(l -^)2 



(28) 



From the root of equation (28), we obtain the first relation between a and 7 as 70 + Oo = 
— ^ — ^y^l + 4: ^^^2j^^^ ■ In a similar fashion, we can obtain other S values and consequently 
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establish a relationship between a„ and 7„, n = 1, 2, 3, • • • as 



X2iz) S2iz) 

Xi{z) si{z) 

Xs{z) S3{z) 

X2{z) S2{z) 



71 + "1 

72 + a2 



3 
2 

5 
2 



2 V 



1 



1 + 4 



D2iii + 1) 
. . . etc. 



(29) 



The nth term of the above arithmetic progression is found to be 

2n + 1 1 



On +7n 



2 V 



(30) 



By substituting for a and 7, we obtain a more explicit expression for the eigenvalues energy as 



nt 



h^Doiii + 1) 



2//i?2 



8/ia2 



yi + 4^^2i!W) + 2^) ' + (?f^ - ^(Z^i + 1)2) 



1^1 + 4^^-^^ + 21^ 



(31) 



Let us now turn to the calculation of the normalized wave functions. By comparing equation 



( 26 ) with equation ( 13 ) we have the following: 

27 + 1 , 2a - 1 , , 

t = ^^, 6 = 1, A^ = -l, m = ^— , a = 2a + l, p = 2(q - 7) + 1. (32) 

Having determined these parameters, we can easily find the wave functions as 

r(2a + n + 1 



-irco 



2F1 (-n, 2(a - 7) + 1 + n; 2a + 1; z) 



(33) 



r(2a + l) 

where F and 2-^1 are the Gamma function and hypergeometric function respectively. By using 



equations (24) and (33), the total radial wave function can be written as follows: 
) + g-iexp(^)^ 



where N^i is the normalization constant 



2F1 -n,2(a-7) + l + n;2a + l; 1 + gexp 



r-R 



(34) 
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Figure 1: The variation of energy spectrum (31 1 as a function of the deformation parameter. For example we 
select II = lfm~^, Vo — 0.3fm~^, R = 7/m and a = 0.65/m. The radial quantum number is fixed to n = 




o 



/= 1 



2 = 2 



Figure 2: The variation of energy spectrum (31 1 as a function of the deformation parameter. For example we 



select n — Ifm ^, Vo — 0.3 fm ^, R — 7/m and a — 0.65 fm. The radial quantum number is fixed to n = 1. 



4 Results and Conclusion 

In Figures ^ and ^ we plot the energy levels Eni versus deformation constant q for ground n = 
and first exited n = 1 states, respectively. In Figure ([T]), it is seen that when q w 0.6, the orbital 
quantum numbers i = 0,1,2, will have the same energy eigenvalues, i.e., E ~ — 0.165/m~^. 
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For the case when q < 0.6, we noticed that Eqq > Eqi > Eq2, whereas when q > 0.6, then 
Eq2 > Eqi > £'00 (less negative or attractive). In Figure Q, taking same parameter set, we 
noticed that when q ~ 1.0, the three states have energy; E -1.04/m-i. In the case when 
q < 1.0, £^10 > £11 > £12 but when q > 1.0, £12 > £11 > £10- The system becomes weakly 
attractive. We have seen that the energy states are sensitive to the deformation constant q. 
There are no physical energy states for I = 0, 1, 2 when g ~ 0.6 for n = and when g < 1.0 for 
n = 1 
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Figure 3: The variation of energy spectrum (311 as a function of the diffuseness of the nuclear surface. For 
example we select n = lfm~^, Vo ~ R = 7/m and q = 1.5 fm. The radial quantum number is fixed to 

n = 0. 



In Figures ([3]) and (|4]), we plot the energy levels E^e versus the diffuseness of nuclear surface 
a for n = and n = 1 states, respectively. For n = case, the energy eigenvalues increasing as 
a increases, i.e., becoming weakly attractive (£02 > Eqi > Eqq). However, for n = 1, the case is 
same as former (i.e. n = 1) but the three curves overlap each other. 

In Figures ^ and (|6j), we plot the energy levels E^e versus the potential depth potential 
depth Vo for n = and n = 1 states, respectively. For n = case, we found that there is a 
substantial change in the energy eigenvalues £"02 > Eqi > Eqq, they increase with increasing Vo 
and becoming less attractive. However, the same behaviour is seen with slow change in three 
curves when n = 1. 

In Figures and ^ we plot the energy levels Ene versus the potential width R for the 
two radial quantum numbers. At i? w 5.1/m, the energy is same (coincides). For R < 5.1/m, 



0.4 



0.3 



4 J 

0.6 




/ = o 



1 



Figure 4: The variation of energy spectrum (311 as a function of the diffuseness of the nuclear surface. For 
example we select jj. = lfm~^, Vo ~ 0.3/m~^, R = 7fm and q = 1.5 fm. The radial quantum number is fixed to 
n = 1. 
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Figure 5: The variation of energy spectrum (31 1 as a function of the potential depth. For example we select 



/i — Ifm ^ , a — 0.65/m, R — 7fm and q = 1.5/m. The radial quantum number is fixed to n = 



Eqo > Eqi > Eq2 (less attractive) but when R > 5.1/m, E02 > Eqi > Eqq. It is obvious that the 
energy levels are sensitive to potential width R. There are no physical energies when R < 5.1 fm 
for n = 0, and R < 6.2/m for n = 1 . For E ^ 0, the particles are no longer attracted. However, 
for £'01 and Eqq —0.2. When n = 1, R ^ 6.2/m, £'02 = £01 = -E'oo- For different values of R, 
the behavior is seen as former for n = 0. 
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Figure 6: The variation of energy spectrum (31 1 as a function of the potential depth. For example we select 
fj, — lfm~^, a — 0.65/m ,R — 0.7 fm and q = 1.5/m. The radial quantum number is fixed to n = 1. 
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Figure 7: The variation of energy spectrum (31 1 as a function of the potential width. For example we select 



/i = 1/m ^ , Vo = 0.3 fm ^, a = 0.65/m and q = 1.5. The radial quantum number is fixed to n = 



In Figures ^ and (10) we plot the energy levels Ene versus the reduced mass /x for the 
two radial quantum numbers. In the ground level the energy becomes repulsive (positive) as n 
increases. For example, E02 becomes positive for /i > Eqi > for fi > and 

Eqo > for /i > 3.0/m~^. This means that the mass has limit. When mass increases, then we 
have no attractive particles and hence no energy spectrum. However, in the excited level, the 
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Figure 8: The variation of energy spectrum (31 1 as a function of the potential width. For example we select 
fj, — lfm~^, Vo = 0.3/m~^, a = 0.65/m and q = 1.5. The radial quantum number is fixed to n = 1. 
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Figure 9: The variation of energy spectrum (31 1 as a function of the particle mass. For example we select 
Vo = O.Sfm~^ , R = 7/m, a — 0.65/m and q — 1.5. The radial quantum number is fixed to n = 0. 



energies of the three states are close to each other and the system remains attractive when the 
reduced mass increases. 

We have seen that the approximately analytical bound states solutions of the wave Schrodinger 
equation for the nuclear deformed Woods-Saxon potential can be solved by proper approximation 
to the centrifugal term within the framework of the AIM. By using the AIM, closed analytical 
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Figure 10: The variation of energy spectrum (31 1 as a function of the particle mass. For example we select 



Vo = 0.3/m ^ , R = 7fm, a — 0.65/m and q = 1.5. The radial quantum number is fixed to n = 1 

forms for the energy eigenvalues are obtained and the corresponding wave functions have been 
presented in terms of hypergeometric functions. 

The method presented in this paper is an elegant and powerful technique. If there are analyt- 
ically solvable potentials, it provides the closed forms for the eigenvalues and the corresponding 
eigenfunctions. However, the case if the solution is not available, the eigenvalues are obtained by 
using an iterative approach [23l [Ml l25] . 
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